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ABSTRACT 

Three epochs of high resolution spectra of the star BD+20 307 show that it is a short period (~3.4 
day) spectroscopic binary of two nearly identical stars. Surprisingly, the two stars, though differing 
in effective temperature by only ^250 K and having a mass ratio of 0.91, show very different Li line 
equivalent widths. A Li 6707A line is only detected from the primary star, and it is weak. This star 
is therefore likely to be older than 1 Gyr. If so, the large amount of hot circumbinary dust must be 
from a very large and recent, but very late evolutionarily, collision of planetesimals. 
Subject headings: stars: individual (BD+20 307) — binaries: spectroscopic — circumstellar matter 



1. INTRODUCTION 

Debris disks older than ~10 Myr containing 
dust at temperatures >100 K are extreme ly rare 
(jAumann fc Probstlll99ll : ISilverstone et al.ll2006[ ). When 
warm dust does appear, it is likely to be from a stochas- 
tic event, perhaps akin to our own solar system's "Late 
Heavy Bombardment," at about 600 Myr after forma- 
tion. Late episodes of dust production may signal the 
presence of a pl anetary system und ergoing architectural 
reconfiguration (jGomes et al.ll2005[) . 

BD+20 307 is one of the few exampl es of a non- young 
star with hot debris (SongetaL 2005) . It has a ring of 
dust at --0.5 AU (Wei nberger e t al. 2008*). In order to 
understand the implications of the large amount of close- 
in dust, it would h elp greatly to know the age of the star. 
I Song et all ()2OO50 used the Li 6707A equivalent width 
and chromospheric activity to suggest an age of ^300 
Myr. 

New observations reported here show that BD +20 307 
is actually a spectroscopic binary. I reexamine the evi- 
dence for the age of the star. 

2. OBSERVATIONS 

The Magellan Inamori Kyocera Echelle (MIKE) spec- 
trograph on the Clay (Magellan II) telescope was used to 
observe BD +20 307 on three consecutive nights - 2007 
October 24-26 (UT). The 0.35" wide x 5" long slit pro- 
vided a resolution of about 55,000 at wavelengths 3400 
- 7250 A. Seeing was '^0.5" on the first two nights and 
0.8-1.2" on the third. On all three nights, data were ob- 
tained with an iodine cell in place to facilitate looking 
for planets around the star, and on the first night an ob- 
servation without the iodine cell was also obtained. I do 
not use the iodine lines for the radial velocity analyses 
that follow. An observing log is given in Table [TJ 

The spectra were flattened, extracted and wavelength 
calibrated using the MIKE pipeline written by D . Kelson 
with methods d escribe d in iKelson et al.l (|2000f ): iKelsonI 
()2003[ ): lKels"oir et al. (2Q0l). The two observations from 
the first night were averaged. The signal-to-noise ratio 
(S/N) per pixel was >100 for wavelengths >4000 A on 
the first two nights, except in the region of maximal io- 
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TABLE 1 
Observing Log 



Date (UT) 


Time 


Int. Time (s) 


2004 Aug 24 
2007 Oct 24 
2007 Oct 25 
2007 Oct 26 


08:28:20 
04:35:20 
04:45:43 
04:12:33 


700 
600 
600 
600 



dine absorption around 5000A. The S/N was about 50% 
worse on the third night due to the worse seeing. 

BD +20 307 was also observed on 2004 August 24 with 
the echelle spectrograph on the 2.5 m du Pont Telescope 
at Las Campanas Observatory. These data cover wave- 
lengths ^4000-9000A and have a resolution of about 
25,000 and S/N of 30-100. The data were extracted and 
calibrated using standard IRAF tasks. 

Heliocentric and barycentric velocities were calculated 
with the RVSAO package in IRAF. 

3. RESULTS 

Two sets of lines are clearly visible in all three nights 
of MIKE data. To obtain the velocities of the double- 
lined spectroscopic binary, cross-correlations with a syn- 
thetic spectrum with effective temperature 6000 K and 
log(g)=5.0 were performed. This spectrum was gener- 
ated using R. O. Gray's SPECTRUM code and line lisfi 
and a Castelli-Kurucz model atmosphere with solar met- 
alicity^. The xcsao package in IRAF was used to com- 
pute the cross-correlations, and the two peaks were fit 
with parabolas in IDL. The uncertainty in the velocities 
was computed as the standard deviation of the velocities 
in the 40 different orders used. Results are reported in 
Table [2] and shown in Figure [TJ On all three nights, the 
primary star produced a higher cross correlation peak. 

The same cross-correlation was done for the lower 
resolution du Pont spectrum. A double-peaked cross- 
correlation appears only for the lowest two orders (4000- 
4150 A). I do not consider this detection of the binary 
reliable, and do not include these RV data in the analy- 
ses which follow. 

^ |http: //www.ph ys . appstate . edu/spectrum/detail .html ' 

2 http://kurucz.harvard.edu/grids/gridP00aODFNEW/ap00k0odfnew.dat 
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TABLE 2 

Measured radial velocities 



Component 2454397.6969 



HJD 

2454398.7041 



2454399.6810 



Primary 
Secondary 



-1.63 ± 0.46 -51.34 ± 0.62 4.92 ± 0.79 
-17.49 ± 0.85 38.03 ± 0.64 -21.24 ± 0.67 





Fig. 2. — A Lithium 6707A line is detected from the primary 
star in all four epochs (solid arrows). These spectra have been 
continuum normalized and offset for clarity. The best upper limit 
(14 mA, 3(t) on the secondary's Li line can be obtained on 25 Oct; 
the dashed arrow shows where it should be located. 



Fig. 1. — Velocities for the primary (filled circles) and secondary 
(open circles) components of the binary phased with an orbital 
period of 3.448 d. The residuals between the fit orbit and the data 
are <1 km s~^ on each point. 



TABLE 3 

Equivalent width (in mA) of Li 6707A line 

Combined Continuum Separate Continuum 
Date Primary Secondary Primary Secondary 



2004 Aug 24 


41 ± 3 


2007 Oct 24 


35 ± 2 


2007 Oct 25 


34 ± 2 


2007 Oct 26 


33 ± 2 



<6 (3o-) 



70 ± 5 
60 ± 4 
58 ± 4 
56 ± 4 



<14 (3cr) 



On all nights a Li 6707A line was detected from the pri- 
mary star (Figure [21) ■ The equivalent widths were com- 
puted using direct integration over the lines relative to 
the combined continua from the two stars. Uncertainties 
from the pipeline reduction were used to give the statis- 
tical uncertainty. An additional systematic uncertainty 
was estimated by choosing different methods of finding 
the continuum and recomputing the equivalent widths. 
A 3(7 upper limit on the secondary's Li line was placed us- 
ing the data from 2007 October 25, when the two stars 
were separated by 89 km a~^ (2 A). These equivalent 
widths are given in the first two columns of Table [31 

The continuum normalized spectrum on 2007 October 
25, which has the maximum separation of the two stars, 
was fit with a combination of synthetic spectra calcu- 
lated with SPECTRUM from Castelli-Kurucz model at- 
mospheres with solar metalicity. Free parameters were 
two effective temperatures, a single wsim, a single log(g), 
and two normalizations. One blue and one red region of 
the spectrum were fit - 4099-4360 A and 6282-6549 A. 
The best fit in both cases had Toff=6500 K and 6250 K 
for the primary and secondary stars, respectively and 
log(g)=5.0. Contours of chi-squarc indicate the uncer- 
tainty is within 250 K (the gridding of the models) in 
Teff. The lines are measurably broader than the ThAr 




3925 



3930 3935 

Wavelength (A) 



3940 



Fig. 3. — Closeup of the calcium K line from 2007 October 25 
showing a weak central reversal for each star. The secondary star's 
reversal (red-ward line) is stronger. The spectrum was continuum 
normalized. 



calibration lamp lines at the same wavelengths. The best 
fit models had vsmi = 10 ± 1 km s~^. 

To compute the stars' Li equivalent widths relative to 
their own stellar continua, the flux ratio of the two stars 
must be obtained at 6707A. Synthetic spectra were fit as 
above to the region at 6645 - 6835A. The best fit flux 
ratio was 1.397 ± 0.007. However, this statistical un- 
certainty probably underestimates the systematics from 
how the models are calculated. Applying this fiux ratio 
to the measured equivalent widths makes the primary 
and secondary equivalent widths increase by factors of 
1.709 and 2.387, respectively. The computed equivalent 
widths are given in the last two columns of Table [H 

The spectra were examined for evidence of chromo- 
spheric activity. Both stars show weak central reversals 
on their Ca H and K lines (Figure [HI that change in ve- 
locity along with the stars. The Balmer lines show no 
central reversals. 
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The velocities were fit for tlie stellar mass ratio (sec- 
ondary/primary) yielding 0.91 ± 0.02 and the line of 
sight radial velocity yielding 7 = —8.8 ± 0.6 km s^^. 

Although the small number of observations prohibit 
the calculation of the full binary orbit, I tested whether 
a circular orbit could fit the velocities. The best fit pro- 
duced reasonable residuals (<1 km s^^) for an orbit with 
a period of 3.448 days, 7=-8.4 km s~^, consistent with 
the fit above, and mass ratio 0.93, again consistent with 
the fit above. This orbit, shown in Figure [1] yields msin"^i 
for each star ~ 0.13 and 0.12 for the primary and sec- 
ondary stars respectively. I estimate the true masses of 
the stars as 1.3 and 1.2 Mq based on the effective tem- 
peratures (6500 and 6250 K) obtained in the spectral fit 
describe d in § 3 and the MK sp ectral-type calibration 
given in lDrilling fc Landoltl ()2000r ): these masses are un- 
certain by ~0.1 Mq. Calculating the inclination from 
these masses yield 27.8 ± 0.8°. 

The fit, assuming the inclination of 28°, gives semi- 
major axes for the binary orbits of 4.4 and 4.7 x 10^ 
km. The stars should have radii of 1.3-1.6 R©, so their 
separation is about 9 stellar radii. 

4. DISCUSSION 

Of primary interest for understanding the dusty debris 
disk around BD-l-20 307 is the age of the system. 

Although the orbit is not well measured, a reason- 
able fit is obtained with e=0. Binaries of any age with 
periods less tha n ^8 days are f ound in circular orbits 
i|Mathieu et"aIlll992l : [Melo et al.ll2001l ). So, e=0 is con- 
sistent with the 3.4 day period and is no indication of 
the age of the system. 

The formation of close binaries may depend on the 
presen ce of another star in the system; iTokovinin et al.l 
(|2006[ ) find that 80% of binaries with periods <7 days are 
actually in higher order multiple systems. It is unknown 
whether BD+20 307 has another, wider, component, but 
this should be investigated. 

Taking sini — 28° from the spectral modeling, the ro- 
tational velocity of the stars is 21 km s~^. A standard 
radius for a Te//=6250 K star would be 1.3 R©, but 
solar-like stars with active chromospheres have some- 
what larger radii than their non-active counterparts 
(jTorreset al.ll2006l ). Stellar radn of 1.3 - 1.6 R© cor- 
respond to rotational periods of 3 - 3.7 d. The close 
agreement between the rotation period and the orbital 
period suggests that the stars are tidally synchronized. 

The timescale for synchro nization for st ars of mass sim- 
ilar to the Sun is given by IZahnI ()1992f ) , and for equal 
mass components is ~ lO^'P^'yr where P is the orbital 
period. Given a period of 3.5 d, tgynch is 100 Myr. This 
provides a lower limit to the age of the system. 

Tidal coupling in a clos e binary can give r ise to en- 
hanced magnetic activity (jRottler et al.ll20d^ . So, the 
central reversals observed in both stars' Ca H and K lines 
need not arise from youth but could be from their inter- 
action. 

For stars with deep convective zones, there is a tight 
correlation between lithium abundance with age and 
spectral type, but for stars with Tg// ~ 6500 K, 
there is little Li deple tion by the age of the Hyades 
(jSoderblom et al.1ll990f ). A typical Li equivalent width 
for a star with B-V of 0.5 in the Hyades is 90 mA, 
1.5 times that for BD+20 307. Furthermore, tidally 



locke d systems in the Hyades show inhibit ed Li deple- 
tion (jBarrado v Navascues fc StauffeHll996l ). 

The situation for older stars is more complex. Stars 
even hotter than the Sun app arently do deplete Li after 
reaching the main sequence. iSestito fc Randich" (20051) 
summarize the Li abundances in open clusters of var- 
ious ages, including the wide spread in abundances 
among members of old (>1 Gyr) clusters. Hot stars 
can show substantial depletions. For example in the 
1 Gyr old cluster NGC 3960, stars with B - V ^ 0.5 
Tef f ~ 6250 K) show a rang e EW(Li) of <10 to 80 mA 
Prisinzano fc Randichll2007n . 

Initial metalicity can have a large effect on main se- 
quence Li abundance, but given their tight orbit, the 
two stars in BD+20 307 presumably formed from the 
same molecu l ar clo ud material at the same metalicity. 
iMartin et all (|2002l ) found that ^20% of wide field bina- 
ries showed disparate Li. Stars may also be polluted by 
Li and other e lements during planetary evolution (e.g. 
lGonzale3 12006[) ■ which could account for differences in 
Li abundance between wide binaries and within clusters 
that should have the same initial metalicity. Especially 
for stars with small convective zones, like the components 
of BD+20 307, material deposited should last for the en- 
tire main sequence lifetime. The tightness of the BD+20 
307 orbit means that any planets must be circumbinary. 
If disparate pollution were to occur, it might be expected 
to favor the more massive star. It is indeed the primary 
star that shows the greater Li line strength. 

Taken together, the modest Li line strength of the pri- 
mary and very low Li line strength for the secondary ar- 
gue for an old age of the system, certainly greater than 
that of the Hyades and probably > 1 Gyr. In this case, 
the dust generating event happened quite late. In the 
Solar System, the last era of major planetesimal removal 
happened at an age 600-800 Myr (Strom et al. 2005). 

It is further tempting to speculate on what happens 
to the dust ring over time. Poynting-Robertson and 
solar wind drag should cause the grains at 0.55 AU 
to fa ll onto th e central stars on timcscales <1000 yr 
CWei nberger et al. 2008) . Meteorites in the Solar System 
contain Li, no matter where they forme d in the original 
proto-planetary disk (jSeitz et al.ll2007[ ). Thus the op- 
portunity seems to exist for substantial stellar pollution 
with Li by the dust around BD+20 307. Perhaps this 
dust is being swept up by a body orbiting between the 
ring and the binary and thus cannot pollute the surface 
of the stars. A detailed abundance study of the stars 
could reveal the extent of any pollution. 

5. CONCLUSIONS 

BD+20 307 is a spectroscopic binary composed of 
two similar stars of late F spectral- type in a short, ro- 
tationally synchronized orbit of ^3.4 days. Standard 
chromospheric-age calibrations for single solar-type stars 
are therefore not applicable to this system. The binary is 
likely quite old, > 1 Gyr, in order to allow the two stars 
to develop quite different lithium equivalent widths. The 
collision that produced the circumbinary dust ring must 
have happened recently, given the short timescale for re- 
moval of material at 0.5 AU, and could be expected to 
pollute the envelopes of the stars. That they have such 
low lithium could be evidence that the dust does not 
actually reach them. 
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